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Abstract 

> 

ff^ ' An analysis of single-electron orbits in combined coaxial wiggler and axial 

o ' 

guide magnetic fields is presented. Solutions of the equations of motion are 

' developed in a form convenient for computing orbital velocity components and 

q ' trajectories in the radially dependent wiggler. Simple analytical solutions are 

• i— i ' 

£>Y obtained in the radially-uniform-wiggler approximation and a formula for the 

43 ' 

derivative of the axial velocity v\\ with respect to Lorentz factor 7 is derived. 
^ Results of numerical computations are presented and the characteristics of the 

equilibrium orbits are discussed. The third spatial harmonic of the coaxial 
wiggler field gives rise to group III orbits which are characterized by a strong 
negative mass regime. 
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I. INTRODUCTION 



Most free-electron lasers employ a wiggler with either a helically symmetric magnetic 
field generated by bifilar current windings or a linearly symmetric magnetic field generated 
by alternating stacks of permanent magnets. A uniform static guide magnetic field is also 
frequently employed. Single-particle orbits in these helical and planar fields combined with 
an axial guide field have been analyzed in detail and have played a role in the development of 
free-electron lasers JIJ . Harmonics of gyroresonance for off-axis electrons caused by the radial 
variation of the magnetic field of a helical wiggler is found by Chu and Lin . Recently the 
feasibility of using a coaxial wiggler in a free-electron laser has been investigated. Freund 
et al. studied the performance of a coaxial hybrid iron wiggler consisting of a central 
rod and a coaxial ring of alternating ferrite and dielectric spacers inserted in a uniform static 
axial magnetic field. McDermott et al. U proposed the use of a wiggler consisting of a coaxial 
periodic permanent magnet and transmission line. Coaxial devices offer the possibility of 
generating higher power than conventional free-electron lasers and with a reduction in the 
beam energy required to generate radiation of a given wavelength. 

In the present paper, single-particle orbits in a coaxial wiggler are studied. The wiggler 
magnetic field is radially dependent with the fundamental plus the third spatial harmonic 
component and a uniform static axial magnetic field present. In Sec. II the scalar equations 
of motion are introduced and reduced to a form which is correct to first order in the wiggler 
field. In Sec. Ill solutions of the equations of motion are developed in a form suitable for 
computing the electron orbital velocity and trajectory in the radially dependent magnetic 
field of a coaxial wiggler. The special case of a radially independent wiggler is also analyzed. 
In Sec. IV the results of numerical computations of the wiggler field components, velocity 
components, radial excursions, and the <3> function for locating negative mass regimes are 
presented and discussed. In Sec. V some conclusions are presented. 
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II. EQUATIONS OF MOTION 



Electron motions in a static magnetic field B may be determined by solution of the vector 
equation of motion 

dv — e „ 

— = v x B (1) 

at '-fine 

where v, — e, and m are the velocity, charge, and (rest) mass, respectively, of the electron. 
Lorentz factor 7 is a constant given by 

7 = (W/c 2 )- 1/2 (2) 

where v = |v| is the constant electron speed. 

The total magnetic field inside a coaxial wiggler will be taken to be of the form 

B = B r .r + B z z, (3) 
B r = B w F r (r,z), (4) 
B z = B + B w F z (r,z), (5) 

where B is a uniform static axial guide field, and F r and F z are known functions of cylindrical 
coordinates r and z. Equation (1) may be written in the scalar form 



^-^ = -^0 + ^), (6) 
at r 

ve + veVr = + q wF ^ _ Vz £} w F r , (7) 
at r 

—j^ = v e Q w F r ; (8) 
Q and Q w are relativistic cyclotron frequencies given by 

n = ^ 0) 

7mc 

V w = (10) 

7mc 

Initial conditions will be chosen such that the transverse motion of the electron in the 
Bq field vanishes in the limit as B w approaches zero. Then, in order to develop a solution to 
first order in the wiggler field B w , the scalar equations of motion will be approximated by 



= tt v r - v\\Vt w F r) (12) 

£ = * < 13 > 

with the wiggler field approximated by the fundamental plus the third spatial harmonic 
component, 

F r = F r isin(k w z) + F r3 sin(3k w z) , (14) 

where 

F rn = G~ l [S n h{nk w r) + ^^(n/c^r)], (15) 
G n = I (nk w R out )K (nk w R in ) - I (nk w R in )K (nk w R out ) , (16) 

2 Ti7T 

= — sin(—)[K (nk w R in ) + K (nk w R out )}, (17) 
nir 2 

2 7Z7T 

T n = — sin(—)[I (nk w R in ) + I (nk w R out )}, (18) 
mr l 

and n — 1, 3; i?j n and _R out are the inner and outer radii of the coaxial waveguide, k w = 2n/\ w 
where \ w is the wiggler (spatial) period, and I , Ii, K , and K x are modified Bessel functions. 



III. ORBITAL ANALYSIS 
A. Radially dependent wiggler 

The scalar equations of motion may be solved to determine the electron orbital velocity 
and trajectory in a coaxial wiggler. Equation (13) yields 

v z = v\\ (19) 

where the constant v\\ is the root-mean-square axial velocity component. With the initial 
axial position taken to be zq = 0, 

z = v\\t. (20) 
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Equations (11), (12), (14), and (20) may be combined to obtain 

cPv, 



dt 2 

where 



+ n 2 v r = fit) (2i) 



fit) = Q Q w v\\ [F r isin(k w v\\t) + F r3 sin(3k w v\\t)]. (22) 

By the method of variation of parameters, a solution of Eq. (21) may be obtained in the 
form 

rt 

v r = [—veo + fi 1 / f{r)cos(n T) drhiniflot) 
Jo 
ft 

+[v r0 -Qn 1 f(r)sin(Q r) dr}cos(tt t) (23) 
Jo 

where v r0 and vq are initial radial and azimuthal velocity components. Then Eq. (11) yields 

ve = No - fi 1 / f( T ) cos (^oT) dT]cos(Q t) 
Jo 

+ [v r0 -fl 1 f(T)sin(floT) dT]sin(fl t) . (24) 
Jo 

The orbital velocity is given to first order in B w by Eqs. (23), (24), and (19). The trajectory 
may then be computed using 

r = r + f v r (r)dT, (25) 
Jo 

= 00+1' v e (r)dr, (26) 
Jo 

and Eq. (20). 

B. Radially uniform wiggler 

By neglecting the radial variation of F r i and F r3 , a solution of Eq. (21) may be obtained 
in the form 

v r = aisin(k w v\\t) + a3sin(3k w v \\t), (27) 

where 



n2_J k 2 v 2 (n = l,3). (28) 

Equation (11) then yields 

vq = — Qq 1 k w v\\aiCos{k w v\\t) — Qq 1 (3k w v\\a 3 )cos(3k w v\\t) . (29) 
The corresponding initial conditions are 

v r0 = 0, (30) 

v eo = -Qq 1 k w v\\a 1 - Q 1 (3k w v \\a 3 ). (31) 

Root-mean-square values of the velocity components may be determined by use of Eqs. 
(27), (28), and (19). Replacing v 2 by its root-mean-square value in Eq. (2) then yields 

V j[l + lQ 2 + \% 2 kla\ + \Q-f + \nfklal] = 1 - 7 - 2 . (32) 

The derivative of v\\ with respect to 7 may be obtained from Eq. (32) and, after some 
algebra, cast into the form 

dv\\ c 2 
dl 77y- 

where 



-$ (33) 



Sn=i, 3 (gg ~ n'klvir^jnlF^nl + 3n 2 klv\) 

2 + E„=i, 3 (ng - n^ w |)-3n2,i^ng(ng + sn^^f ) ■ 1 j 

This equation may be used to establish the existence of a negative mass regime. 



IV. NUMERICAL RESULTS 

A numerical computation is conducted to investigate the properties of the equilibrium 
orbits of electrons inside a coaxial wiggler. Wiggler wavelength 2ir/ k w and lab-frame electron 
density no were taken to be 3 cm and 10 12 cm~ 3 , respectively. The wiggler magnetic field B w 
was taken to be 3745 G which corresponds to the relativistic wiggler frequency Q w /ck w = 
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0.442. Electron-beam energy (7— 1)itlqc 2 was taken to be 700 keV, corresponding to a Lorentz 
factor 7 = 2.37. The axial magnetic field B was varied from to 25.3 kG, corresponding to 
a variation from to 3 in the normalized relativistic cyclotron frequency Qo/ck w associated 
with Bq. The inner and outer radii of the coaxial wiggler were assumed to be Ri n = 1.5 cm 
and R^t = 3 cm, respectively. 

Figure 1 shows the variation of the axial velocity of the quasi-steady-state orbits with 
the axial guide magnetic field for three classes of solutions. Group I orbits for which < 
flo < k w v\\, group II orbits with k w v\\ < Q$ < 3k w v\\, and group III orbits with Q > 
3k w v\\. Existence of group III orbits is due to the presence of the third spatial harmonic of 
the wiggler field, which also produces the second magnetoresonance at Qq ~ 3k w v». The 
narrow width of the second resonance at Qo/ck w m 2.7 compared to the width of the first 
magnetoresonance at fi ~ k w vu is illustrated in Fig. IB. This is due to the relatively weak 
third harmonic compared to the fundamental component of the wiggler field. It should be 
noted that although the exact resonances f2 = k w vu and Qq = 3k w v n occur at the origin 
where Vn/c— Qo/ck w = 0, the first "magnetoresonance" in the literature refers to the group 
II orbits with the cyclotron frequencies around VL^/ckw ~ 1 in Fig. 1. Similarly we refer 
to the group III orbits with the cyclotron frequency around fl /ck w ~ 2.7 as the second 
magnetoresonance . 

The rate of change of the electron axial velocity with electron energy is proportional to $ 
and is equal to unity in the absence of the wiggler field. Figure 2 illustrates the dependence 
of $ on the radial wiggler magnetic field and the axial guide magnetic field Bq. The curves 
corresponding to the group I and II orbits are almost unaffected by the third harmonic and 
are almost the same as in Ref. where the third harmonic is neglected. A negative mass 
regime (i.e., negative $ for which a decrease in the axial velocity results in an increase in 
the electron energy) is found for group III orbits which is stronger than that of the group II 
orbits. 

Equations (15)-(18) are used to calculate the radial components of the wiggler field F rl 
and F r 3. For the axial component the following expressions are used [|| 
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F z = F zl cos(k w z) + F z3 cos(3k w z), 



(35) 



F zn = G n l [S n I Q (nk w r) - T n K (nk w r)]. (36) 

Figure 3 shows the variation of the amplitudes of the wiggler magnetic field (divided by 
B w = 3745 G) with radius, for the first and third spatial harmonics. For the first harmonic 
the radial component, F r \, has a minimum at r « 2.28 cm and the axial component, F z i, 
changes sign around this point. Mc Dermott et al. || have demonstrated the stability of 
a thin annular electron beam when F r \ is minimum at the beam radius. The radial and 
axial components of the third harmonic of the wiggler F r3 and F z3 are also shown in Fig. 3. 
Magnitudes of F r3 and F z3 both are minimum at r ~ 2.28 cm where F r \ is minimum. This 
is actually an inflection point for F z3 . 

Variations of the radial components of the first and third spatial harmonics of the nor- 
malized wiggler magnetic field F r ± and F r3 with the wiggler wave number k w are shown in 
Fig. 4. Figure 4 also shows the dimensionless transverse velocity coefficients a\ = a\jc and 
a 3 = a 3 /c for the initial orbit radius r ~ 2.28 cm where F r \ is minimum. The cyclotron 
frequency VL /ck w « 2.7 is taken at the second magnetoresonance, and our choice of 3 cm 
for the wiggler wavelength corresponds to k w « 2.1 cm" 1 . It can be observed that at this 
wave number although the radial component of the wiggler field at the first harmonic F r x 
is much larger than the third harmonic F r3 , the transverse velocity coefficients of the third 
harmonic a 3 are larger than oc\. This shows that the third harmonic may have considerable 
effects around the second magnetoresonance at Qq ~ 3k w v». Away from this resonance Eq. 
(28) shows that a 3 will be of the order of F r3 . 

In order to study the transverse motion of electrons in the radially dependent wiggler 
field Eqs. (11), (12), (25), and (26) are solved numerically with the initial conditions chosen 
so that, in the limit of zero wiggler field, there is axial motion at constant velocity vu but 
no Larmor motion. Figure 5 shows the variation of the radial and azimuthal components of 
electron velocity with z (= vut). The normalized cyclotron frequency Qo/ck w is chosen to be 
0.5, 1.2, and 3 for group I, II, and III orbits, respectively, which are somewhat away from the 
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magnetoresonances. Solid curves correspond to the initial orbit radius r = 2.28 cm, which 
is at the point where F rl is minimum. Broken curves correspond to r = 1.8 cm, which is 
away from the F r \ minimum. It can be observed that the spatial periodicity of v r and v$ for 
the first two groups is equal to one wiggler wavelength, which is the same as that of the first 
harmonic. Although group III orbits have a clear sinusoidal shape at r = 2.28 cm (solid 
curves), the slight deviations from sinusoidal shape are obvious at r = 1.8 cm (broken 
curves). This is because that at r = 2.28 cm where F ri is minimum F r3 is very small. 
Therefore away from the resonance the third harmonic plays almost no role, at r = 2.28 
cm. Moving away from F r i minimum to tq = 1.8 cm, however, increases the magnitude of 
F r3 slightly making the effect of the third harmonic noticeable on group III orbits, which 
are away from the second magnetoresonance, in Fig. 5. 

Figure 6 shows the variations of v r /c and vg/c with z (= v\\t) for group III orbits when the 
cyclotron frequency is adjusted at the second magnetoresonance at Qo/ck w ps 2.7. At r = 
2.28 cm the periodicity is approximately equal to \ w /3, which is the same as that of the third 
harmonic, and shows the strong influence of the third harmonic on the transverse velocity 
components. Going away from the F rl minimum to r = 1.8 cm makes the amplitudes of 
oscillations of v r and v e larger. Broken curves correspond to the solutions of Eqs. (27) and 
(29) when the radial variation of the wiggler field is neglected. These solutions do not differ 
appreciably from the r-dependent solutions for r = 2.28 cm because at F rl minimum the 
radial excursions are small for group III orbits. Away from the F rl minimum at r = 1.8 
cm, however, deviations are noticeable due to the larger radial excursions. 

Figure 7 shows v r /c versus z for group II orbits for r-dependent wiggler (broken curves) 
and r-independent wiggler (solid curves). The initial orbit radius is taken at r = 2.28 cm 
and the cyclotron frequency is chosen around the first magnetoresonance at fl /ck w = 0.9. 
Large radial excursions of electrons for group II orbits make the transverse velocity strongly 
affected by the radial dependency of the wiggler field. The amplitude is also modulated in 
space with the wavelength of around 16X W = 48 cm. 

The radial excursion r shown in Fig. 8 corresponds to the cyclotron frequencies away 
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from the magnetoresonances at Qo/ck w equal to 0.5, 1.2, and 3 for the group I, II, and III 
orbits, respectively. Solid curves correspond to the r = 2.28 cm and the broken curves 
correspond to To — 1-8 cm. It can be noticed that when the electrons are injected into the 
wiggler at ro = 2.28 cm, where F r % is minimum, electron orbits remain well away from the 
waveguide walls at Ri n = 1.5 cm and R out = 3 cm. 

Figure 9 compares the radial excursions of group III orbits at the second magnetores- 
onance at Qo/ck w = 2.7, (solid curve) with those slightly away from the resonance at 
Qo/ck w = 3, (broken curve). Influence of the third harmonic can be clearly seen through 
the modulation of the third harmonic by the first harmonic when the cyclotron frequency is 
adjusted at the second magnetoresonance. 

V. CONCLUSIONS 

The third spatial harmonic of the coaxial wiggler field gives rise to the group III orbits 
with Q > 3k w v\\. This relatively weak third harmonic makes the width of the second 
magnetoresonance narrow compared to the first magnetoresonance. A strong negative mass 
regime is found for the group III orbits. By adjusting the cyclotron frequency at the second 
magnetoresonance the wiggler induced velocity of the group III orbits was found to increase 
considerably. When the electrons are injected into the wiggler where its magnetic field is 
minimum the electron orbits remain well away from the waveguide boundaries. 

Harmonic gyroresonance of electrons in the combined helical wiggler and axial guide 
magnetic field is reported by Chu and Lin . In their analysis the relativistic single particle 
equation of motion is used with the axial velocity as well as the axial magnetic field of 
the wiggler averaged along the axial direction. By assuming near steady-state orbits for 
off-axis electrons they found that the radial variation of the wiggler magnetic field produces 
a harmonic structure in the transverse force. This force, in turn, comprises oscillations at 
all harmonics of k w z. It should be noted that there is no harmonic structure in the helical 
wiggler itself and the higher velocity harmonics vanish for the exact steady-state orbits of 
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the on-axis electrons. Moreover, higher harmonics do not appear in the one dimensional 
helical wiggler where the radial variation is neglected. 

In the present analysis of coaxial wiggler, on the other hand, equation of motion is 
written to first order in the wiggler amplitude. With this approximation axial component of 
the wiggler field has no contribution to the problem leaving the axial velocity as a constant. 
Magnetic field of a coaxial wiggler is composed of a fundamental plus a large number of odd 
spatial harmonics, which directly appear in the magnetic force represented by f(t) in Eq. 
(22). Third harmonic in f(t) appears in the transverse velocity components as a part of 
the integrands in Eqs. (23) and (24) and is also demonstrated numerically for the radially 
dependent coaxial wiggler, but for the radially uniform wiggler the third harmonic is explicit 
in the solutions Eqs. (27)-(29). 
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FIGURES 

FIG. 1. Normalized axial velocity vn/c as a function of the normalized axial-guide magnetic 
field flo/ck w for group I, II, and III orbits. Narrow width of the second resonance at £lo/ck w 2.7 
compared to the width of the first magnetoresonance at fio ~ k w v\\ is illustrated in Fig. IB. 

FIG. 2. Factor # as a function of the normalized axial-guide magnetic field Qo/ck w for group 
I, II, and III orbits. 

FIG. 3. Radial dependence of the radial and axial magnetic fields (divided by B w = 3745 G) 
in the coaxial wiggler for the fundamental and third spatial harmonics. 

FIG. 4. Wave-number dependence of the radial components of the normalized wiggler mag- 
netic field F r i, F r 3 and the dimensionless transverse velocity coefficients «i, Q3. The normalized 
cyclotron frequency Qo/ck w = 2.7 is taken at the second magnetoresonance and ro = 2.28 cm. 

FIG. 5. Normalized transverse velocity components as a function of axial distance z for the 
initial orbit radius ro = 2.28 cm (solid curves) and ro = 1.8 cm (dotted curves). The normalized 
cyclotron frequency Qo/ck w is 0.5, 1.2, and 3 for the group I, II, and III orbits, respectively. 

FIG. 6. Normalized transverse velocity components as a function of axial distance z for group 
III orbits, at the second magnetoresonance £lo/ck w = 2.7, solid curves correspond to the radial 
dependent wiggler and broken curves correspond to the radial independent wiggler. 

FIG. 7. Normalized radial velocity as a function of axial distance z for group II orbits at the 
first magnetoresonance O,o/ck w = 0.9 with ro = 2.28 cm. Solid curves correspond to the radial 
dependent wiggler and broken curves correspond to the radial independent wiggler. 

FIG. 8. Radial excursion r as a function of axial distance z for ro = 2.28 cm (solid curves) and 
ro = 1.8 cm (broken curves) for groups I, II, and III orbits out of resonance at £lo/ck w =0.5, 1.2, 
and 3, respectively. 

FIG. 9. Radial excursion r as a function of axial distance z for group III orbits at ro = 2.28 
cm. Solid curve corresponds to resonance at Qo/ck w = 2.7 and broken curve corresponds to out of 
resonance at Qo/ck w = 3. 
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